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i , Abstract 

The Bariev model with open boundary conditions is introduced and analysed in 

CJ | detail in the framework of the Quantum Inverse Scattering Method. Two classes 

of independent boundary reflecting i^-matrices leading to four different types of 

boundary fields are obtained by solving the reflection equations. The models are 

exactly solved by means of the algebraic nested Bethe ansatz method and the four 

sets of Bethe ansatz equations as well as their corresponding energy expressions are 

derived. 
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1. Introduction 
O 



" ^^ i 

^ Much work has been devoted, during the last years, towards a better under- 

standing of integrable models of strongly correlated electrons. Such an interest 
arises from the efforts to unveil a theory that can provide a consistent picture 
for observed phenomena of condensed matter physics such as high-T c super- 
conductivity. Two prototypes of integrable models for this purpose are the 
one-dimensional (ID) Hubbard model solved by Lieb and Wu [1] (see also [2]) 
and the t-J model [3] (which is integrable for the super symmetric coupling [4]), 
mainly because of their relevance in the description of electronic mechanisms. 
This is justified on physical grounds since the electron hopping is strongly 
disturbed by the on-site Coulomb interaction in the ID Hubbard model and 
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by the spin fluctuations through the antiferromagnetic coupling for the super- 
symmetric t-J model. More recently, motivated by the inclusion of additional 
interactions whether through internal impurities or external boundary fields, 
a renewed interest in the study of integrable strongly correlated electron sys- 
tems has taken place. Some examples of the latter are the ID Hubbard open 
chain [5], the supersymmetric t-J model with boundary fields [6,7] and other 
integrable electronic systems with boundary interactions associated with Lie 
superalgebra symmetries [8]. 

The study of integrable systems with open boundary conditions (OBC) [5,9- 
11] has turned into an active domain since the introduction by Sklyanin [12] 
of a generalization of the Quantum Inverse Scattering Method (QISM) [13,14] 
which provides a systematic approach to handle the boundary problem by 
introducing a set of equations, called reflection equations (RE). The solutions 
of these equations, referred to as the boundary K-matrices, in turn introduce 
boundary interactions into the Hamiltonian of the system in such a manner 
that integrability is preserved. Indeed, the presence of boundary fields dras- 
tically changes surface critical properties, and in some occasions results in 
the appearance of so-called boundary bound states [15-17]. Nevertheless, in 
spite of these results, several models still lack a comprehensive study regarding 
OBC, one of those being the ID Bariev model [18,19]. 

The ID Bariev (interacting XY) chain is also a Hubbard-like integrable model 
of special interest, as it exhibits the existence of hole pairs of Cooper type 
which are relevant to theories of superconductors. The integrability of the ID 
Bariev model within the framework of the QISM was first established by Zhou 
in [20] (see also [21]). To achieve this goal one has to build non-trivial higher 
conserved currents that allow the identification of a quantum i?-matrix satis- 
fying the Yang-Baxter equation [13,14]. This .R-matrix does not possess the 
difference property for the spectral parameter, a feature that is also shared 
by the i?-matrix associated with the Hubbard model [22]. Soon after this 
parametrization of the model had been obtained, the algebraic Bethe ansatz 
solution for the periodic chain was calculated independently by Martins and 
Ramos [23] and by Zhou [24]. On the other hand, as far as we know, studies 
concerning OBC for the Bariev model appear incomplete. This can be under- 
stood in light of the fact that the i?-matrix is presented in a complicated form 
with the absence of both additivity for the spectral parameter and crossing 
unitarity. As such, we must resort to searching for solutions of a much more 
intricate RE. Though some known solutions to the RE for the ID Bariev open 
chain provide a certain class of boundary K-matrices [25] , a thorough investi- 
gation is still lacking and, moreover, the algebraic Bethe ansatz solution has 
not been obtained yet. 

In this paper, we shall study in detail the OBC for the ID Bariev chain 
within the framework of the QISM. We will show that besides the .ftf-matrices 



obtained in [25] there exists a new class of solutions of the RE which leads to 
pure magnetic boundary fields in the Hamiltonian and provides more general 
nesting boundary A'-matrices for a hidden asymmetric XX Z open chain. Here 
it is worth noticing that this class of boundary fields may have a feasible 
realization by applying boundary external voltages in experiments on quantum 
wires [26]. By means of the algebraic Bethe ansatz technique, we obtain the 
eigenvalues of the transfer matrix of the ID Bariev model for four different 
kinds of possible boundary fields. For all cases we present the Bethe ansatz 
equations and the energy spectrum in explicit form, which provides a first step 
towards investigating the surface critical properties for the ID Bariev model 
with these different boundary fields. 

The paper is organized as follows. In section 2 we construct the Bariev model 
with four kinds of boundary fields by means of the QISM adapted to special 
boundary conditions. The basic quantities, e.g., the .R-matrix, the matrices K± 
defining the boundaries, the monodromy matrices and the transfer matrices 
are defined. In particular, we solve the RE and obtain two independent classes 
of solutions for the boundary A'-matrices using the variable-separation method 
[27]. The relation between the transfer matrices and the Bariev Hamiltonians 
with different boundary fields is established. Sections 3 and 4 are devoted to 
the solution of the models through Bethe ansatz methods. In particular, in 
section 3 we discuss the action of the transfer matrix on the pseudovacuum 
state. In section 4 we perform the nested algebraic Bethe ansatz and find 
the Bethe ansatz equations, the eigenvalues of the transfer matrices and the 
energy spectrum of the system. Section 5 is reserved for our conclusions and 
discussions. 



2. The Bariev open chain with boundary fields 



Let us begin by introducing the Bariev model with boundary fields whose 
Hamiltonian reads 



JV-l 



H = E {(4^+i + a i a Hi) ex p (/^Vi+i" 

+ (t?t- +1 + t-t? +1 ) exp (#7+07) } + b.t. (2.1) 

where the four different boundary terms are given by 



b .t.(i,i) = _^M [sinh j3a z N r z N + cosh/3 (a z N + r z N )} 
C2 



exp Gg) 



[sveih. (3 o{tI + cosh/3 (a{ + rf)] , 



b.t. (n ' n) = \ exp (/3) c 2 a z N + \ exp (/3) Cl rf , 
b t (i,ii) = _^Pj^ [ginh/?(TM + cosh/? ^ + ^j 

+ -exp( / 3)c 1 r 1 2 , 
b.t. (II ' I} = ^ exp (P) c 2 a 



N 



(>XP ( ^ [sinh/3a z rf + cosh/5 (<r* + rf)] . 



Cl 



(2.2) 
(2.3) 

(2.4) 

(2.5) 



Above o~i and Tj are the two commuting species of Pauli matrices acting on 
site i, f3 is a coupling constant and C\, c 2 are the boundary parameters char- 
acterizing the strength of the boundary fields. Notice that in cases (2.2,2.3) 
the boundary terms at the left and right ends are consistent with the bulk 
symmetry which is a combination of the inversion j — > N — j + 1 and the 
exchange a ^ r. The cases (2.4,2.5) are interchanged under this symmetry. 

In order to derive these models, we first consider the following i?-matrix [20] 

R(u, v) = 
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(2.6) 



whose elements (Boltzmann weights) are given explicitly in Appendix A. 
It satisfies the non-additive Yang-Baxter equation 

R 12 (u, v)R 13 (u, w)R 23 (v, w) = R 23 (v, w)R 13 (u, w)R 12 (u, v) (2.7) 



and enjoys the following properties 



R 12 (u,v)R 2 i(v,u) 
Rl\(-v,u)R%(u,-v) 

R%(-u,v)Rl\(v,-u) 



1, 

1. 



(2.8) 

(2.9) 

(2.10) 



Above, tj denotes matrix transposition in the ith space and R 2 \ = P.R\ 2 .P, 
where P is the permutation operator. The matrix R is a new object introduced 
in Ref. [25] in order to handle integrable boundary terms. We may take either 
(2.9) or (2.10) as its definition. We remark here that neither the crossing 
unitarity nor the difference property is valid for this i?-matrix, which turns 
the calculations into a complex problem. 

The Yang-Baxter algebra can also be verified for the monodromy matrix T 



Ru{u, v) T (u) T (v) =T (v) T (u)R 12 (u, v) 



(2.11) 



which is defined, as usual, by 

T(u) = L m (u) ■■■L 01 (u). 



(2.12) 



Above the subscripts and N, ... ,2,1 denote the auxiliary and quantum 
spaces, respectively and the local Lax operator is given by 



Lj(u) 



/i(u)<t;qj 

\h(u)f 2 (u)a+ 



h{u)P+rr f 1 (u)a~Q+ h{u)f 2 {u)ajT-\ 

P+Qj f!(u)a-r+ f 2 (u)a-Qj 

ft{u)a+T- PfQ+ f2(u)Pfrr- 

r+ f2(u)a+Qj f 2 (u)p-T+ P 3 Q~ J 



(2.13) 



where 



n 



pf 



Q 



2 V 3 ' 

-(l±a 
2 

1 



ii 



-exp(/3)(l T (Tj); 



j- 



-(i=F<#; 

exp(/3)(l T r/); 

ITT/); 
h{u) = v /l + W 2 exp(2/3), f 2 (u) = v / TT 






W 



Next, following the QISM adapted to the case of special boundary conditions, 
we define the doubled monodromy matrix as 



U-(u) = T{u)K_{u)T~ 



■u 



(2.14) 



such that the transfer matrix is given by 

T (u) = Tr K+(u)UJu) 



(2.15) 



Above T 1 is the inverse of the monodromy T and K± are the matrices defining 
the boundaries. The requirement that they obey the RE [25] 



R 12 (u,v) K_ (u)R 21 (v, -u) K_ (v) 

2 1 

—K- (v)Rn(u, -v) K_ (u)R 2 i(-v, — u), 
Rl\ t2 (v,u) K% (u)R 12 (-u,v) K^ (v) = 
K^ (v)R 21 (-v,u) K% (u)R\\ t2 (-u,-v) 



(2.16) 



(2.17) 



together with the Yang-Baxter algebra assure that the transfer matrix com- 
mutes for different spectral parameters, proving the integrability of the model. 
Therefore the transfer matrix (2.15) may be considered as the generating fuc- 
tion of infinitely many integrals of motion for the system. We emphasize that 
there is no isomorphism between the matrices K + (u) and K_(u), due to the 
absence of the crossing unitarity for the i?-matrix. Therefore, we have to solve 
the two RE separately in order to fix the boundaries. In Appendix A, this 
calculation is presented in detail using the variable-separation prescription. 
This approach, which was already proposed in [27] to treat the D 2 model 
with OBC, is extended here for the case of an .R-matrix without the difference 
property. It is in fact a systematic method which simplifies drastically the 
calculation of the K± matrices. After a tedious algebra, we find two differ- 
ent classes of boundary K± matrices consistent with the integrability of the 
Hamiltonian (2.1), listed below 



K + (u) 



(K\±\ 




V o 



u 




K2±l 
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KZ±\ 
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(2.18) 



KA ± (u)) 



Solutions for K_, corresponding to the left boundary 
Case I 



Kl_{u) = (u — ci){h 2 u — Ci), 



K2^(u) = (ci - u)(h 2 u + ci), 

K3_(M) = (c 1 -M)(/ i 2 M + Ci), (2.19) 

K4_(w) = (u + Ci)(/l 2 M + Ci), 

Case II 

/a_(u) = (l + ci/m), 

K2_(m) = (1-ci/im), 

A"3_(u) = (l + ci/i«), (2.20) 

K4_( y u) = (l-c 1 hu). 

Solutions for K + , corresponding to the right boundary 
Case I 

Kl + (u) = (c 2 u + l){h 2 c 2 u + 1), 

K2 + {u) = h 2 {c 2 u-l){c 2 u + l), 

K3 + (u) = (c 2 u-l)(c 2 u + l), (2.21) 

KA + (u) = (c 2 u - l)(c 2 u - h 2 ), 

Case II 

K\ + {u) = h(hu — c 2 ), 

K2 + (u) = h(hu-c 2 ), 

K3+(u) = (u + hc 2 ), (2.22) 

K4 + (u) = (u + hc 2 ) 

where h = exp {(3). These two classes of boundary if-t-matrices provide four 
possible choices of BC, according to the combination of the boundary pairs 
(K + (u),K-(u)), which are being labeled by (1,1), (II, II), (I, II) and (II, I) and 
originate the boundary terms b.t. , b.t. ' ', b.t. ' and b.t. ' ' respec- 
tively. Taking into account the fact that TrK + (0) = and TrK,(0) = in 
the cases (1,1) and (I, II) and TrK + (0) = in the cases (II, II) and (II, I), we 
can show that the Hamiltonian (2.1) with the boundary terms b.t/ ' ' and 
b.t.' ' ' is related to the transfer matrix matrix (2.15) in the following way 

r(u) = C x u 2 + C 2 (H + const.) w 3 + • • • , (2.23) 

while 

T (u) = C 3 u + C 4 (H + const.) u 2 + • • • (2.24) 

for the Hamiltonian (2.1) with the boundary terms b.t/ ' ' and b.t/ ' '. The 



difference in the leading term expansions (2.23,2.24) is a consequence of the 
fact that for case I the matrix K + (u) is quadratic in the variable u while 
for case II it is linear. Above Cj, % — 1, • • -4, are some scalar functions of the 
boundary parameters. In the next sections, we shall be focusing in the solution 
of the eigenvalue problem of the transfer matrix (2.15) and, consequently, in 
the diagonalization of the present model. 



3. Algebraic Bethe ansatz I: Action of the operators on the pseu- 

dovacuum 



In order to carry out the algebraic Bethe ansatz for the Bariev model with 
boundaries, we first need to find the eigenvalues and eigenvectors of the trans- 
fer matrix (2.15) 

r|$„) = A|$ n ) (3.25) 



This will be achieved by extending the Bethe ansatz techniques developed in 
[28], [29] for the Hubbard model to the present case. Following this prescription, 
the eigenvectors of the transfer matrix are given by 

|$ n ) = * n .F|0). (3.26) 



where the components of F are coefficients of an arbitrary linear combination 
of vectors <fr n (see next section for details) and |0) is the pseudovacuum state, 
choosen here as the standard ferromagnetic one 

|0) = ®f =1 |0)„ (3.27) 



where 



|0>i= ' i ® " (3-28) 



which corresponds to the doubly occupied state. Then it is easy to work out the 
action of the monodromy matrices T, T _1 on this state. We begin by writing 




the periodic monodromy matrix T(u) as 



T{u) 



B{u) B 1 (u) B 2 (u) F{u) 

Ci(w) A n (u) A 12 (u) E l {u) 

C 2 (u) A 21 (u) A 22 (u) E 2 (u] 

\C 3 (u) C 4 (u) C 5 (u) D{u) ) 



(3.29) 



and T(u) = T 1 (—u) as 



f{u) 



( 



u 



\ 



u 



B{u) B 1 (u) B 2 {u) F[ 

Ci(«) A n (u) A 12 (u) E v 

C 2 (u) A 21 (u) A 22 (u) E 2 (u) 

\c 3 (u) C 4 (u) C 5 (u) D(u) ) 



(3.30) 



By computing the action of T(u) and T(u) on the psedudovacuum we find 



[«rio>, 



lN 



B(«)|0> = B(«)|0> = |0>, 
D(u)\0) = D(u)\0) 

A aa ( W )|0) = A aa (u)\0) - [u exp (P)} iy |0), 
A 12 («)|0)=^ 12 («)|0) = 0, 
A 21 («)|0)=^ 21 («)|0) = J 
S o (u)|0)^0, B a (u)\0)^0, 
E a (u)\0)^0, E a (u)\0)^0, 
F(u)\0)^0, F(u)\0)^0, 

c t (u)\o) = a t (u)\o) = o, 

i — 1, • • • , 5, a, b — 1,2. 



(3.31) 



Similarly, the doubled monodromy matrix LL = T(m)X_(w)T(m)(2.14) can be 
arranged as 



tf_(«) 



/ B(u) Si(u) B 2 («) F(u) \ 

c^w) in(u) ii 2 («) ^i(«; 

c 2 (n) l 21 («) l 22 («) £ 2 ( w ; 



(3.32) 



9 



Then using eq.(3.31) and the Yang-Baxter algebra 

2 2 

T (u)R 12 (u, -u) T (u) =T (u)R 12 (u, -u) f (u), (3.33) 

we can calculate the action of the doubled monodromy [/_ on the pseudovac- 
uum, which due to the existence of two classes of boundary K_ matrices (2.19), 
(2.20) yields (see Appendix B) 

B(u)\0) = Wr(u)B(u)B(u)\0), (3.34) 

A aa (u)\0) = W~ +1 (u)A aa (u)A aa (u)\0), (3.35) 

D(u)\0) = W 4 (u)D(u)D(u)\0}, (3.36) 

B o («)|0)^0, E a (u)\0)^0, 

A ab (u)\0)=0, F(u)^0, 

Ci(u)|0) = 0, i = l,---,5, a^b= 1,2, (3.37) 

where for case I: 

wr(u) = i, 

, N . , 2h 2 u(uci - 1) 

w 2 („) = w 3 (») = - (ftV _ 1)(tf „_ Cir 
w , = 2 M 2 ( UCl -i)(; i 2 + i)( 1Kl -; i 2 ) 

(w 2 — 1) (-u 2 — h 2 ) (h 2 u 2 — c\) (u — C\) ' 
and for case II: 

wr(u) = i, 

ur-t ^ 2/lV 



2r; (/i% 2 -l)' 

_ 2/m(/m + Ci) 
3 W ~ (/i% 2 - 1) (1 + c x hu) ' 

2 M 3 (fe 2 + l)( M + fe Cl ) 
4W (« 2 -l)( M 2 -/i 2 )(l + Cl /m)' 

Actually, in order to considerably simplify the calculations, the following trans- 
formations 



A aa {u) = A aa {u)- p2 f> U ( b(u 
Pi{u,-u) 

D(w) = D(-u) t \B(u) 

pi(u,-u) 



' 1 
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1+U \ln(u)-^±^A 22 (u), (3.38) 



h? — u 2 h 2 — u 

were used in the above equations (3.34)-(3.36). In fact, these types of transfor- 
mations have also been performed in other models with OBC [5,9-11]. They 
simplify not only these expressions for the action of the operators on the pseu- 
dovacuum, but also simplify the fundamental commutation relations, allowing 
for an easy recognition of the "wanted" and "unwanted" terms, essential to 
the Bethe ansatz machinery (see next section for details). Then, after this 
transformation and taking into account the two sets of boundary matrices K + 
(2. 21), (2. 22), the transfer matrix (2.15) can be recast as 

t{u) = W+{u)B{u) + J^ W+ +1 (u)A aa (u) + W+(u)D(u), (3.39) 

o=l 

where for case I: 

2u 2 (u - c 2 ) (h 2 + 1) {h 2 u - c 2 ) 



wr(u) 



W+(u) 



W 3 + (u) 



Wt(u) 



(uc 2 + 1) (uh 2 c 2 + 1) (u 2 - 1) (h 2 u 2 - 1) 

u (uc 2 - 1) (u - c 2 ) (h 2 + 1) 
(h 2 uc 2 + l)(u 2 -h 2 )(uc 2 + \y 
h 2 u (uc 2 - 1) (u - c 2 ) (h 2 + 1) 
(h 2 uc 2 + l){u 2 -h 2 )(uc 2 + iy 
(uc 2 - 1) (uc 2 - h 2 ) 
{uc 2 + I) (h 2 uc 2 + l)' 



and for case II: 

Wf(u) 
W+(u) 

WUu 



2u (h 2 + 1) (hc 2 u + 1) 
h{u 2 -l){h 2 u 2 -l){hu-c 2 ) 
{h 2 + l)(u + hc 2 ) 
h {u 2 — h 2 ) (hu — c 2 ) 
u (h 2 + 1) (h + uc 2 ) 



' ' {u 2 - h 2 ) {hu - c 2 ) 

(u + hc 2 ) 



W£{u 



h {hu — c 2 ) 



Notice from eqs. (3.34)-(3.37) that the Lax operators (3.32) act as creation 
(annihilation) operators on the pseudovacuum. Before ending this section, we 
would like to recall that the boundary matrix pairs (K + (u), K_{u)) lead to 
four kinds of boundary fields preserving the integrability of the model. Con- 
sequently, the factors W^{k) i = 1, ... ,4 also form a pair (W^~(k), W~{k)) 
leading to four possible choices consistent with the boundary matrix pairs. The 
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perfect factorization provide us with the nesting boundary X-matices appear- 
ing in (3.35) and (3.39) respectively. In the following section, we shall prove 
that they indeed constitute the solutions of the RE for the hidden asymmetric 
XX Z spin open chain. 



4. Algebraic Bethe ansatz II: The nesting procedure 



The next step in the solution of the eigenvalue problem of the transfer matrix 
consists in the computation of the necessary commutation relations between 
the diagonal fields and the creation fields. From the RE (2.16) and the defini- 
tion (3.32), after tedious calculations, we obtain the main fundamental com- 
mutation relations. We display below just those needed for future calculations 



B(u)B a (v) = 
D{u)B a {v) = 
A adl (u)B d2 (v 



(u + v)(uvh 2 + 1) ~ 



(u — v)(uvh 2 — 1) 
(u+ v)(uvh 2 + 1) 5 



B a (v)B(u)+u.t., 
B a {v)D(u)+u.t., 



(u — vh 2 )(uv — 1) 
(u + v)(l + uvh 2 ) 
(u — v){uvn z — 1) 



d\d 2 
b 2 



xB C2 (v)A Clbl (u) + u.t., 

B a (u) ® B b (v) = {B c (v) <g> B d (u) 



+- : : F(v)i [I <g> A(u)) \ .r(-v, -u) 



1 — uv 

Vi + h 2 v 2 VTT 



u z ~ 



1 — uv 



(h 2 u 


+ 


v) y/\ + h 2 U 


2 Vl + V 2 r 


1 


(u 

+ 


— h?v) {uv - 


-1) [J 


(h\ 


u) y/\ + h 2 V 


Vi + « 2 5 



{h 2 u + v) Vi + h 2 u 2 VTT 



F(u)Z[I®A(v) 

F(v)B(u) 

F(u)B(v) 



■*, 



(4.40) 

(4.41) 

(4.42) 



(4.43) 



where 



£ = (0,/r\l,0), A(u) 



A u (u) A 12 {u) 
A 21 (u) A 22 {u) 



(4.44) 



12 



r 12 (u,-v) 



r(—v, —u) 



i 



^10 0^ 

c(u, —v) d(u, —v) 
a(u, —v) c(u, —v) 
1 / 

10 0^ 

a(—v, —u) c(—v, —u) 
c(—v, —u) d(—v, —u) 
\0 ly 



(4.45) 



(4.46) 



and 



a(u, —v) = 
d(u, —v) = 
c(u, —v) = 
a(—v, —u) 
d(—v, —u) 
c(—v, —u) 



pr{u, -v)pi(u, -v) - p$(u, -v) 

pl(u, -v) - pl(u, -v) 
p 10 (u, -v)pi{u, -v) - f$(u, -v) 

p\{u,-v) - (%(u,-v) 
Pi(u,-v)p 15 (u,-v) 



pl(u,-v) - pl(u,-v) 

= Pu(-v, -u)p 10 (-v, -u) - pj 3 (-v, -u) 

Pi(-v,-u)p 14 (-v,-u) 
= Pu(-v, -u)p 7 (-v, -u) - pl 2 (-v, -u) 

Pi{-v,-u)pu(-v,-u) 
__ Pu{-v , -u)pi 5 (-v, -u) - pu(-v, -u)p 13 (-v, -u) 
Pi(-v,-u)pu(-v,-u) 



. (4.47) 



It remains to consider the commutation relations for the creation field F. 
Since they are extremely involved, we have collected them in Appendix B. We 
remark here that the fundamental commutation relations are much more com- 
plicated in comparsion with those obtained in the case of periodic boundary 
conditions, turning the calculations here much more cumbersome and diffi- 
cult to treat. In all relations, we have omitted the unwanted terms (those 
that can not generate an eigenvector of the transfer matrix) to save space. 
It turns out that the auxiliary matrices f 12(11, —v) and r(— v, — u) are noth- 
ing but the rational R matrices of an isotropic asymmetric six-vertex model. 
Indeed, by performing the parameterization u = exp (iki) ,v = exp (ifc 2 ) and 
h = exp (/?) = exp (— irj), one may present these matrices as 
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r(—v, —u) = r(k\ — k 2 ) 



^10 

a(ki — k 2 ) c(ki — k 2 ) 

c{h - k 2 ) d{ki - k 2 ) 

1 
f 12 (u, -v) = P.r(k 1 + k 2 - 27/), 



V 



r 21 (u, -v) = P.r(ki - k 2 ).P 
with the Boltzmann weights 



(4.48) 



(4.49) 
(4.50) 



a(k) 



d(k) 



c(k) 



sin 



sin rj ( 1 



2 i '/ 

suit/ 



sm 



(| + ,) eXP (^)' 



sin f ( ■ 1 u\ 



sm 1 1 + 77 



(4.51) 
(4.52) 

(4.53) 



In view of the commutation relation (4.43), we may phenomenologically write 
[28] the n-particle state as 



\$ n (k 1 ,...,k n ))=$ n (k 1 ,...,k n )F\0). 



(4.54) 



with a recursive relation of n-particle vector 



$ n (A4,..., & n ) = B ai (ki) ®$ n -l(k2,...,k r , 



it 

+ E { [t® ^( fc l)] ®n-2(k 2j ..., kj-x, kj+i 



i=2 



xB{k 3 )Gf\{k u ...,k n )} 

n 
+ E { [I® ^( fc l)] ^n- 2 (fc 2 , • • • , fcj-1, %+l 

i=2 

x [/®i(A;,)]i/j™ ) 1 (A;i,...,A; n )}. 



, . . . , fv n J 



, . . . , /c n j 



(4.55) 



Above F = F ai ' • • • ' fln are the coefficients of arbitrary linear combination of the 
vectors reflecting the 'spin' degrees of freedom with <2j = 1,2 (i = 1, . . . ,n) 
and £ plays the role of forbidding two spin up or two spin down particles at 
the same site. F(u) creates a local particle pair with opposite spins. From the 
commutation relation (4.43), we can also conclude that $ n (&i, • • • , k n ) satisfies 
the symmetry relation 
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^n\Ki, " ' " , fcj, Kj+1, ' ' ' ; Kn) — 

<&„(fci, • • • , kj+i, kj, ■■■ , k n ).r(kj - k j+1 ) (4.56) 

This symmetry giving a restriction to the functions iJJ_ 1 (/ci, • • • , k n ) and 

GjZ x (k-\_, ■ ■ ■ , fc n ) is very useful to deduce the coefficients and simplify the un- 
wanted terms in the eigenvalue of the transfer matrix. In fact, after performing 
three-particle scattering, the explicit form of these coefficients can be fixed. 
However checking three-particle scattering is indeed a extremely tough calcu- 
lation and we had to leave the coefficients to be determined later. Explicitly, 
we display the two-particle state 

$ 2 (h,k 2 ) = B ai (h) <g> B a2 (k 2 ) 
J I + e i2fcl \/l + e i2( - k2 ~^ -> 

+ [_% 1+k2) — nh) i (i®A( k2 )) 



cos(^^ + r])Vl + e i2k Wl + e i2 ^-^ ~ 
s i n (fci_M + 7? )(l_ e i2(fci+te)) 



F(h)B(k 2 ). £ . (4.57) 



Implementing a rescaling k — * k + r\ and seting c\ = e l2 ^~ , c 2 = e l2 ^+ , the 
eigenvalue of the transfer matrix (2.15) acting on the state (4.54) can be given 
as 



TL pac I "'~rfej 



r(fc)|$ n (fci, . . . , M) = W+(k)Wr(k) I] 



cos I — j* ) COS 



(¥0 



i=l 



sin f^ 1 ) sin ( 
+ JlWt +1 {k)W- +1 {k)^ Nk 

a=l 

n cos f m^) cos (^0 , , 

X TT MrT TT^^aa (k, {h}, {AO,-}) 

t X i sin (*±k) sin (*=**) aa V ' X "' 1 3S) 

n C oS ( fc+fc ' ] COS ( k ~ ki ] 

+W+(k)W^(k)e [m( ^ k+ ^ J] ' V 2 ; ' V 2 ^ 



2 



i=i sin (^ + rfj s i n (i_k + ^ 
|$„(A:i,...,A;„)) (4.58) 

provided that 

^(SS^) '* - »= - A( " (* = *■■ W- M) ■ (4-M) 

i — 1,- ■ -n. 

Here A^ (A;, {/cj}, {Aj}J is the eigenvalue of the reduced transfer matrix r^(k), 
which turns out to be an inhomogenous isotropic six-vertex model with open 
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boundary conditions, namely, 

r (1 \k)F ai --- a " =A il) {k,{k l },{X J })F ai --- a ", (4.60) 

where 

T (1) (fc) = Tr K { +\k)T^(k)K { 2\k)T^'\-k). (4.61) 

The nesting monodromy matrices T^ l \k) and T^ (— k) read 

TW(fc) = r 12 (fc + fci)£V 12 (* + fc)^ • • -r 12 (fc + fc,,)^. (4.62) 
TW _1 (-fc) = r 21 (fc - /c n fe^" • • .r 21 (fc - fc)££r 21 (* - h)f ev (4.63) 

and the nesting K±(k) may be chosen as 
Case I 

K?(k)=\ . I, (4.64) 



e 



ii) 



in such a way that the factors W + can be rewritten as 



W + (k) = sin(| + §-e + )sin(|-f-e + )oos>y 

11 J cos(| + f + ^ + )cos(|-§+e+)sinA:sin(A: + ?7)' l ' J 

W + (k) - W + (k) - Sin( ^ + I + e+) sin( * + i ~ e+) COST? f4 66) 

M/ 2 (*) - M/ 3 (k) - ms{k _ _ v _ + ^ +) cQs{k _ + v _ + ^ +) s . n(A; + ^ , (4.66) 



sin(| + f + e + )sin(| + f + ^ 
cos(f + § + £ + )cos(f -§+£+) 



Wt{k) = ""VLTVT Tl + , (4-67) 



and 



^'(fc)= | 10 | (1.68) 

1 



with 



Wf(fc) = 1, (4.69) 
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wv(*)= .,» T^^T^:^,, „, (4.7i; 



sm(| - I - £_) sin(| + f - f_) sin(ft + 77) sin(ft + 2r/) : 



Case II 



'W/ 



cos(| + 77 - £ + )e 



-ik 



K?(k)= ^ |, (4.72) 

\ C0S(| +7/ + £ 



with 



and 



, . , cos(| + £+) cosr^e lfe 

^ fc = ■ Tfe \ ; . , ^77- v (4-73) 

sm(| — £ + ) sin ft sm(ft + 77) 

W+(ft) = W+(ft) = - . (h ^77—^ v (4-74) 

sm(| - 4+) sin (ft + 277) 

, , s cos(| + 77- £+)e ir ? 
Wt ft = l2 ' .7 , 4.75 



m , cos § + £_ e lfc 

Ki 1} (ft) = I V2 ^ I , (4.76) 

cos(§-£_ ' 



with 



^f(ft) = 1, (4.77) 

W 2 (k) = W 3 ~(k) = . - ) k , (4.78) 

sin ft cos (2 + ?-) 

W7(ft) = cos (f+^-^) cos ^ i(fc+??) (4.79) 

cos(| + £_) sin(A; + 77) sin(A; + 277) 

We notice that the nesting boundary K± '( ft) -matrices also constitute four 
classes of BC consistent with the boundary matrix pairs (K + (k), K~(k)). So 
far, the eigenvalue problem of the ID Bariev model with four kinds of pos- 
sible boundary fields reduces to the solution of the nested auxiliary trans- 
fer matrix (4.60), which can be associated with an inhomogeneous isotropic 
six- vertex model with open boundary conditions (4.64)-(4.76). Furthermore, 
we can check that the r-matrix (4.48) and the nesting K± -matrices (4.64)- 
(4.76), which arises from the factorization of the spin degree of freedom, indeed 
satisfy the Yang-Baxter algebra 

12 21 

r 12 (ft - A) T« (ft) T« (A) =T<U (A) T« (ft)r 12 (ft - A), (4.80) 
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and the reflection equations 



1 2 

-(1) <U\„..<U , U ITU) 



r 12 (/c - A) Ki i; (k)r 21 (k + A) Ki 1J (A) 



2 1 



^^ (A)ri 2 (fe + A) K y y (k)r 21 (k - A), (4.81) 



1 2 



r^ t2 (A - fc) ^ x (fc)ri 2 (-fe - A) i^ i; (A) 

2 1 

=K^ t2 (A)f 21 (-fc - A) K? tx (A:)r t 1 1 2 t2 (A - fc), (4.82) 



where 



r 12 (k)r 21 (-k) = l, (4.83) 

^ 2 i(-fcM a 2(fc) = l, (4-84) 

r\l(-ky 2 \(k) = l. (4.85) 

Following all steps in solving XXZ open chain in [12], one can completely 
diagonalize the nesting transfer matrix (4.61) with the help of the main ingre- 
dients (4.80) -(4.82) describing the open BC for the hidden asymmetric XXZ 
open chain. Therefore, what follows is the exact solution of the XXZ open 
chain characterizing the spin degree of freedom, i.e., 

Al'>(fc,{fc,},{A,})|* (1) (A 1 ,...,A„)) 

sm{k + rj) j=l sin (_j. + n) sm(-^) 

sm{k + i]) j=1 sm (__j--|- 7 yj sm (_j.j 

x |$«(A^--,A M )>, (4.86) 

provided that 

WMk)Wn(k)A^(k) " s jn(^ + 27y)sin(^ + T y) 

W^(fc)W^(*)I>«(*) U sin(^)sin(^-r/) ' l ' j 

1 = 1, 

where 

AW(Jfc) | $ (1) ) =| $ (1) ), (4 



sm('*±*Osin(*=^ 



D u(k) i $«) = n ,r+? 2 »v x i $(1) >- ( 4 - 89 ) 

VJl ' \\\ sin(^i + 77) sin(^i + 77) J ' 7 l ^ 

Here for case I, we have 



W+(fc) = W^k) = 1, (4.90) 

W+(fc) = W^(fc) = 1, (4.91) 

while for case II, 



W+(k) = cos(^ - £ + )e~ ik , W u (k) = cos(^ + ^)e ifc , (4.92) 

W+(k) = cos(^ + r) + £+)<*, W£(fc) = cos(^ + 77 - £_)e-\ (4.93) 

respectively. 

Finally, making a shift on the spin rapidity Aj — »■ Aj — 77, the eigenvalue of the 
transfer matrix (2.15) with two classes of boundary K-matrices can be given 
explicitly (up to a common factor) 



n cosf^ 1 ) cosf^-^' 



r(«) 1 * n (h, . . . , k n)) = |n sm (A; sin (A^ {wmwr(k) 

W?{k) W^{k) W& (k) W n (k) e i: 



^ n ( k + 2 V) Jxr +n.\jjr-n.\j Xr +n.\jjr-n.\„i2Nk 



sin(/c + 77) 

M 



M sin(^^2) sin(^^ 



2 / v 2 
sml^^smf 



1=1 °"H 2 J S1I H 2 / 



( . ()S /fe±fci') C0S (fc_Ji' 



t t=i sin(^ + r/) sin(Y + ?7) L 

SU) *' -H/3 + (A;)H/3-(fc)M/ 2 + 2 (A;)H/ 2 - 2 (fc)e i2JVfc 



sin (A; + 77) 
'n ^;^,' ))!^,^)), (4.94, 

where all factorized coefficients Wf~ (i — 1, • • • , 4, 11 and 22) have been given 
explicitly in the above equations, which provide four kinds of possible exact 
solutions. Finally, the expression above is in fact the eigenvalue of the transfer 
matrix, if the following Bethe ansatz equations are satisfied for the charge and 
spin rapidities ki, Xj . 
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e-%(^ + )C 2 (^-)^n : in i4^:J44 (4-95) 

sin( J 2 ' — 77) sin( J 2 ' — 77) 
t = 1, 

j = l,...,M, i = l,...,n. 

The expressions for Ca(&z>£±) an d ^(Aj,^), a = 1,2 depend on which type 
of boundaries we are considering. Below we display all cases 
Case (1,1) 



rin( ! p+e+ U(fc,^)=^- ( f t 

sin(^2-e+) sin(V-^ 



&(*.,*+) = "TV" r , fc(fc,,«-) = ""A^- , (4-97) 



;^2-e+)' sin(^p 

^ 1 (A„e+)=^(A J ,^) = l, (4.98) 

Case (II,II) 



ci(m+) = - cos( i' H » Ca( ^'^ } = - gifci ' (499) 

cos(f + £+) 

^(A J ,€,) = C " ( ^-y . (4.100) 

cos(^+£ + ) 

^■t-)= -% + y . (4.10D 

COS ( 2 ~ f+) 



Case (1,11) 



sin(^±a+^ 
sin(^2 - £ H 



Ci(fe,£+) = rU% 7; . C 2 (^,C-) = "A (4.102) 



coGf AjT; + E-)e lXi 

M*j,t+) = 1, V*(A*£-) = 1+, ; ; , (4-103) 

cos(^ 4+) 



Case (11,1) 



cos(» - g+) ... p ,_sin(^ + e-; 
• cos( ^ + e+ )'^'^--( 



ci(*.,w = -^Ttv^v &<*.*-) = :,% t , ( 4104 ) 
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M*s,t + ) = ^ % ]+ri c : , M^Z-) = 1. (4-105) 

cos(^+£ + ) 

Above M is the number of electrons with spin down and n is the total number 
of the electrons. Considering in addition the relations (2.23) and (2.24) one 
can find the energy spectrum for all cases 
Case (1,1) 

E = - C -^- C -^ + 4±cosk t , (4.106) 



C\ c 2 



i=\ 



Case (II,II) 

n 

^ = exp(/3)(c 1 + c 2 )+4^cosA; i , (4.107) 



i=l 



Case (1,11) 



cosh 8 n 

£ = exp(/3) Cl -^^ + 4^008^, (4.10? 

C 2 i=l 



Case (11,1) 



cosh 8 n 

E=-^^- + exp(p)c 2 + 4Ycosk i . (4.109) 



Conclusion and discussion 



In summary, we have analysed in detail the open boundary conditions for the 
Bariev chain with special boundary conditions. Two classes of boundary K + - 
matrices have been obtained by solving the RE. It has been found that these 
two classes of solutions of the RE lead to four kinds of integrable boundary 
fields for the charge and spin degrees of freedom separately. Through the 
nesting procedure, we have diagonalized exactly the two level transfer matrices 
with four kinds of possible boundary fields. The eigenvalues of the transfer 
matrices, the energy spectrum and Bethe ansatz equations have been derived 
for all cases. The boundary fields indeed contribute nontrivially to the ground- 
state properties as well as the low-lying spectrum. The functions Ca(^«, £±) and 
■0a(-\/,£±) w iH contribute nontrivial phase-shift factors to the density of the 
charge rapidity and spin rapidity. We notice that the Bethe equations (4.95) 
and (4.96) would be reduced to the purely doubling of the ones for the ID 
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Bariev model with periodic BC [23] [24] if the boundary parameters £± — > oo. 
These boundary parameters £±, denoting the impurity strength too, might 
change the band filling, the boundary surface energy, the mesoscopic effects 
as well. The results obtained provide us with a basis to investigate the surface 
critical properties, and correlation functions [14,30,31] for the resulting model. 
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Appendix A. The Boltzman weights and the boundary .K\|--matrices 



Before practising our ansatz for fixing the boundary K±(u)-m&trix, we display 
the Boltzmann weights of the .R-matrix (2.6), which appeared first in [20] 



p2 ~- 
Pa ~- 
P5 = 
P% ~- 
Ps ~- 
Pio 
P12 

Pl4 



{u — v)h 
1 + uvh 2 

Vi + /» 2 uVi + /iVvTTuVTTu 5 



VI + h 2 u 2 Vl + hh 
1 + uvh 2 



-Pi, P3 



Pi, 



[1 + uv)(l + uvh 2 ) 



-Pi, 



hy/l + h 2 v 2 \/l + u 2 {u - v) 
(1 + w)(l + uvh 2 ) 

Vl + h 2 v 2 Vl + u 2 (u - v) 
(1 + uv)(l + uvh 2 ) 

VTT^VTTv 2 



pi, 



Pl, P7 



h 2 (u — vf 



1 + uv 



-Pi, P% 



(1 + uv)(l + uvh 2 ) 
(u — v)(u — h 2 v) 



Pi, 



[1 + uv)(l + uvh 2 ) 



Pi, 



(u — vY 



(1 + uv)(l + uvh 2 ) 
hy/l + h^VTTv 2 ^ - v) 
(1 + uv)(l + uvh 2 ) 
(u — v){h 2 u — v) 



Pi, Pu 



(u — v)h 
1 + uv 



-Pi, 



[1 + uv)(l + uvh 2 ) 



Pi, Pl5 



Vl + h 2 u 2 y/T+v*{u - v) 
(I + uv)(l + uvh z ) 
h(u — v ) 2 



uv)(l + uvh 2 ) 



Pi, 



As mentioned in section 2, due to the absence of the isomorphism between 
X-t-matrices, we have to solve the RE (2.16) (2.17) separately to fix them. 
Substituting K-(u) (2.18) into the RE (2.16), we find, at first glance, that 
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the RE (2.16) involve two variables u and v which make the functional equa- 
tions involving Ka_(u), a = 1,...,4 much more complicated. However, if 
we observe the structures of the i?-matrix and the RE (2.16), we may pick 
up some simpler functional equations of the RE which allow us to separate 
Ka-(u),Ka-(v) into the following forms 



P2(u, v)p 3 (v, -u)Kl_(u) + p 3 (u, v)p 2 (v, -u)K2_(u) 
p 2 (u, -v)p 3 {-v, -u)Kl_(u) + p 3 (u, -v)p 2 (-v, -u)K2_(u) ' 

p 2 (u,v)p 3 (v,-u)Kl_(u) + p 3 (u,v)p 2 (v,-u)K2_(u) 
p 2 (u, -v)p 3 (-v, -u)Kl_(u) + p 3 (u, -v)p 2 (-v, -u)K2_(u) ' 
pu{u, -v)p 8 (-v, -u)K3-(u) + p 8 (u, -v)pu(-v, -u)KA_(u) 

p u (u, v)p 8 (v, -u)K3_(u) + p 8 (u, v)p 2 (v, -u)KA_(u) 

Substituting the Boltzmann weights of the i?-matrix into above equations, and 
analying the consistant conditions of these equations, we can easily conclude 
an ansatz 



K\_ 


(v) 


K2_ 


(v) 


K\_ 


(v) 


K3_ 


(v) 


K3_ 


(v) 


KA_ 


(v) 



Kl_(u) 

K2„(u) 
K3-(u) 
KAju) 



■■ (ci + c 2 u)(c 3 + c 4 u)(c 5 + Cqu), 
(ci - c 2 u)(c 3 + c A u)(c b + c 6 u), 
(ci + c 2 u)(c 3 - c A u)(c b + Cqu), 
(ci + c 2 u)(c 3 - c 4 u)(c 5 - ceil), 



(A.l) 
(A.2) 
(A.3) 
(A.4) 



with minimal coefficients q to be determined. Running the RE again with this 
ansatz (A.1)-(A.4), it is easily found that only one coefficient is free. Thus the 
two classes of boundary /^.-matrices can be immediately chosen as the forms 
(2.19) and (2.20). 

Simlarly, substituting fT + -matrix (2.18) into the RE (2.17) we have 



Kl + (v) 



p 2 (v, u)p 16 (v, -u)Kl + (u) + p 3 (v, u)p 3 (v, -u)K2 + (u) 



K2 + (v) p 3 {-u, -v)p 2 (u, -v)Kl+(u) + p 2 (-u, -v)p 16 (u, -v)K2 + (u) ' 
K\+(v) p 2 (v,u)p 16 (v,-u)Kl + (u) + p 3 (v,u)p 2 (v,-u)K3+(u) 



K3. 
K3. 



[V) 

(v) 



p 3 (-u, -v)p 3 (u, -v)Kl+(u) + p 2 (-u, -v)p 16 (u, -v)K3 + {u) ' 
p 8 (-u, -v)p 15 (u, -v)K3+(u) + pu(-u, -v)p 7 (u, -v)KA + {u) 



KA + (v) p n (v, u)p 9 (v, -u)K3+(u) + p 8 (v, u)p 15 (v, -u)KA + (u) 

here we first have to determine p's from the relation (2.9) and (2.10), i.e. 



P2(x,y) 



Ps(x,y) 



(1 + h z xy)\\ + xy)y/\ + h 2 yWl + h 2 x 2 


(x — y) 2 (h 2 x — y)(h 2 y — x) 


(1 + h z xy) z (l + xy)y/l + Wy'Wl + h 2 x 2 



(x — y) 2 (h 2 x — y){h 2 y — x)h 2 
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P7(x,y) 

P9(x,y) 
Pi5(x,y) 
Pw(x,y) 



(1 + h l xy)\\ + xyf^l + y 2 - 


sj\ + x 2 


(x — y) 2 (h 2 x — y)(h 2 y — 
(1 + h 2 xy) 2 {\ + xyfh 2 sj\ + • 


x) 
i/Vl + x 2 



(x — y) 2 (h 2 x — y)(h 2 y — x) 
[1 + h 2 xy)(l + xy) 2 



(x — y) 2 (xh 2 — y) 
(1 + h 2 xy) 2 (l + xy) 
(x — y) 2 (yh 2 — x)h 



With the help of the new Boltzmann weights, the consistency of above equa- 
tions indeed can be able to lead to an ansatz 



Kl + (u) 
K2 + (u) 
K3+(u) 
KAju) 



(ci + c 2 u)(c 3 + C4uh 2 )(c 5 + c 6 u), 
(c 2 u - h 2 ci)(c A uh 2 + c 3 )(c 5 + cqu), 
(ci + c 2 u)(c 4 u - c 3 )(c 5 + ceii), 
(ci + c 2 u)(c 4 u - c 3 )(c 6 u - c 5 h 2 ). 



(A.5) 
(A.6) 
(A.7) 
(A.8) 



Running the second RE (2.17) again, the solutions (2.21) and (2.22) would be 
fixed. 



Appendix B. Useful commutation relations 



For factorizing the transfer matrix (2.15) acting on the pseudo- vacuum state, 
we need the following commutation relations 



C 1 (u)B 1 (u) 
C 1 (u)B 1 (u) 
C 2 {u)B 2 {u) 
C 2 (u)B 2 (u) 
Cs(u)F(X) 



^4 \B{u)B{u 
Pi(u,-u) L 

p 2 (u,-u) r - 

— BlujBlu 

Pi{u,-u) L 

P2JU, ~U) r- 

— r B(u)B(u 

pi{u,-u) L 

P2JU, -U) r 

pi(u,-u) 

P4JU, -U) 

Pi(u,-u) 

P2JU,-U) 

Pi(u,-u) 



B(u)B(u 

B(u)B(u 



A u (u)Au(u) 
A 22 (u)A 22 (u) 
A 22 (u)A 22 (u) 
D(u)D(u 



C 3 (u)F(u) = P ^ U \ c 1 (u)B 1 (u) + 

p 2 {u, —u) 



C 4 (u)£i(u) + C 5 (u)E 2 {u 
p 8 (u, -u) 



;b.i) 
;b. 2 ) 
;b.3) 
;b.4) 



p 2 {u, -u) 



A u (u)A u (u) 



(B.5) 
D(u)D(u) 
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pi(u, -u) =, p 7 (u,-u) , =, 

: C 4 («)£?i(«) 7 -C 5 (w)£ 2 (w), (B.6) 



p 2 (u,-u) 
C 3 {u)F{u) = P ^~ U \ c 2 (u)B 2 (u) 



p 2 (u, -u) 
p 10 (u,-u) 
p 2 (u,-u) 



CAujEAu 



p 2 (u,-u) 
Ps(u, -u) 
p 2 (u,-u) 
_ pi(u,-u) 
p 2 (u, -u) 



A 22 {u)A 22 {u) - D{u)D{u 
C 5 (u)E 2 (u), (B.7) 



which can be derived directly from (3.33). In addition, one also can show the 
relations Ci(X)B a (X) = for i ^ j, i = 1,2,3, and a = 1,2. Using these 
commutation relations and after lengthy algebra, one may manage to present 
the factorized forms of (3.35) and (3.39). 

In order to check the eigenvalue problem of two-particle excitation, the fol- 
lowing commutation equations are necessary 



+Pi 2 (-v,-u 
+Pn{-v, -it 
+pu(-v,-u 



Pi(u, v)pu(v, -u)F(u)B(v) = p u (u, -v)p 4 (-v, -u)F(v)B(i 
' ■ ' l p 3 (u, -v)B 2 (v)Bi(u) + p 12 (u, -v)F 2 (v)A u (u 

p 3 (u, -v)Bi(v)B 2 (u) + pn(u, -v)F 2 (v)A 22 (u 

Pl (u, -v)B{v)F{u) + p 2 (u, -v)fli(u)£i(u) 

+p 2 {u, -v)B 2 (v)E 2 {u) + p 4 {u, -v)F(v)D(u)] , 

P9(u,v)p u (v, -u)D(u)F(v) + 

p 5 (u,v) [p 12 {v, -u)i4n(u)F(v) +/9n(u,-u)^i(«)£2(v) 

+p 6 (u,v) [p 13 {v, -u)A 22 (u)F(v) + p u (v, -w)S 2 (w)^i(^) 

+Pa{u,v) p 4 (v,-u)B(u)F(v) + p 8 (v, -^B^ujE^v) 

+p$(v, -u)B 2 (u)E 2 {v) + Pl (v, -u)F(u)D{v) 

= Pa{u, -v)pi(-v, -u)B(v)F{u) + pi(u, -v)p!(-v, -u)F(v)D(u) 

+p s (u,-v) pi(-v, -ujB^E^u) + p!(-v, -u)B 2 (v)E 2 (u) , 

p 3 (u,v) \p 3 {v, -u)A u (u)F(v) + p 13 (v, -m)^i(m)E 2 (v) 
+p 2 (u,v) [p 2 (v,-u)B(u)F(v) + p±(v, -u)!}! (u)^ (v) 
+p w (v, -u)B 2 {u)E 2 {v) + p 8 (v, -^F^Div) 1 



Pu(-v,-u) 



p 13 (u, -v)S 2 (v) j Bi(m) + p u (u, -v)F(v)A n [ 



u 



+p$(-v, -u) [p 2 (u, -v)B(v)F(u) + pi(w, -v )B 1 (v)E 1 
+Pio(u, -v)B 2 (v)E 2 (u) + p 8 (u, -v)F(v)D(u) , 
p 3 (u, v) [p 3 (v, -u)A 22 (u)F(v) + p 12 (v, -ujE^E^v) 
+p 2 (u,v) \p 2 (v,-u)B(u)F(v) + p 7 (v, -u)-Bi(it)£i(i;) 



u 



(B.8) 



(B.9) 



(B.10) 
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+p 1 (v, -u)B 2 {u)E 2 {v) + p 8 (v, -u)F(u)D(v) j 
= Pii(-v, -u) p 13 (u, -v)B 2 (v)Bi(u) + pu{u, -v)F(v)A 22 (u 
+Ps(-v, -u) [p 2 {u, -v)B(v)F(u) + p 7 (u, -v)Bi(v)^i(w) 
+01 fw. -v^Mi^fw) + o«fw. -v)F(v)D(u) . 



+pi(u, -v)B 2 (v)E 2 (u) + p 8 {u, -v)F(v)D{v 
Pz{u,v) p 3 (v,-u)A 12 (u)F(v 
= Pn(-v,-u) 



+pi 2 (v,-w)i±a(w)-Ei(? 

L Pi2(u, -u)Bi(v)S 2 (w) + pn(«, -v)F 

Pi 5 (w, -u)P8(-v, -«)B 2 (v)-Ei(m), 

■'" "> rp 3 ( 1 , i _ M )A 21 ( M )F(v) + p 13 (i;, -u)E 2 {u)E 2 (v) 



(B.H) 
+ p 2 (w, v)pi5(w, -m)-B 2 (m)^i(v) 

»il 2 («)" 



p 3 {u,v 



= Pn(-v,-u) p 12 (u,-v)B 1 (v)B_ 
+pu(u, -v)ps(-v, -u)Bi(v)E 2 (u). 



(B.12) 

j +/52(U,U)P15(V,-W)5i(u)^2(v) 

i(u) + Pii(w,-v)F(u)i 2 i(«) 



(B.13) 
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